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DYNAMICS OF THE TIME TO THE MOST RECENT COMMON 
ANCESTOR IN A LARGE BRANCHING POPULATION 

By Steven N. Evans-*^ and Peter L. Ralph^ 
University of California at Berkeley and University of California at Davis 

If we follow an asexually reproducing population through time, 
then the amount of time that has passed since the most recent com- 
mon ancestor (MRCA) of all current individuals lived will change as 
time progresses. The resulting "MRCA age" process has been stud- 
ied previously when the population has a constant large size and 
evolves via the diffusion limit of standard Wright -Fisher dynamics. 
For any population model, the sample paths of the MRCA age pro- 
cess are made up of periods of linear upward drift with slope 
punctuated by downward jumps. We build other Markov processes 
that have such paths from Poisson point processes on R++ x R++ 
with intensity measures of the form A ® /i where A is Lebesgue mea- 
sure, and n (the "family lifetime measure" ) is an arbitrary, absolutely 
continuous measure satisfying ^((0,oo)) — oo and iJ,{{x,oo)) < oo for 
all a; > 0. Special cases of this construction describe the time evolu- 
tion of the MRCA age in (1 -I- /3)-stable continuous state branching 
processes conditioned on nonextinction — a particular case of which, 
/3 = 1, is Feller's continuous state branching process conditioned on 
nonextinction. As well as the continuous time process, we also con- 
sider the discrete time Markov chain that records the value of the 
continuous process just before and after its successive jumps. We 
find transition probabilities for both the continuous and discrete time 
processes, determine when these processes are transient and recurrent 
and compute stationary distributions when they exist. Moreover, we 
introduce a new family of Markov processes that stands in a relation 
with respect to the general (1 + /3)-stable continuous state branching 
process and its conditioned version that is similar to the one between 
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the family of Bessel-squared diffusions and the unconditioned and 
conditioned FeUer continuous state branching process. 

1. Introduction. Any asexually reproducing population has a unique 
most recent common ancestor, from whom the entire population is de- 
scended. In sexually reproducing species, the same is true for each nonre- 
combining piece of DNA. For instance, our "mitochondrial Eve," from whom 
all modern-day humans inherited their mitochondrial DNA, is estimated to 
have lived around 180,000 years ago [19] while our "Y-chromosomal Adam" 
is estimated to have lived around 50,000 years ago [33]. There have also 
been efforts to estimate the time since the MRCA lived (which we will also 
call the "age of the MRCA") in populations of other organisms, particularly 
pathogens [32, 34]. These studies, using sophisticated models incorporating 
of demographic history, are focused on estimating the age of the MRCA at 
a single point in time (the present). 

As time progresses into the future, eventually the mitochondrial lineages 
of all but one of the daughters of the current mitochondrial Eve will die out, 
at which point the new mitochondrial Eve will have lived somewhat later 
in time. The age of the MRCA is thus a dynamically evolving process that 
exhibits periods of upward linear growth separated by downward jumps. 

Recently, [26] and [29] independently investigated the MRCA age process 
for the diffusion limit of the classical Wright-Fisher model. The Wright- 
Fisher model is perhaps the most commonly used model in population dy- 
namics: each individual in a fixed size population independently gives birth 
to an identically distributed random number of individuals (with finite vari- 
ance), and after the new offspring are produced, some are chosen at random 
to survive so that the total population size remains constant. The diffusion 
limit arises by letting the population size go to infinity and taking the time 
between generations to be proportional to the reciprocal of the population 
size. 

In this paper, we investigate the MRCA age process for a parametric 
family of population models in a setting in which the population size varies 
with time, and, by suitable choice of parameters, allows control over the 
extent to which rare individuals can have large numbers of offspring that 
survive to maturity. The model for the dynamics of the population size is 
based on the critical (1 + /3) -stable continuous state branching process for < 
/3 < 1. These processes arise as scaling limits of Galton- Watson branching 
processes as follows. 

(n) 

Write for the number of individuals alive in a critical continuous 
time Galton- Watson branching process with branching rate A and offspring 
distribution 7. The distribution 7 has mean 1 (and thus, the process is 
"critical"). Suppose that if is a random variable with distribution 7, then 
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the random walk with steps distributed as the random variable {W — 1) falls 
into the domain of attraction of a stable process of index 1 + /3 G (1,2]. The 
case /3 = 1 corresponds to 7 having finite variance and the random walk 
converging to Brownian motion after rescaling. Set xj.^"^ = n~'^I^Z^^ and 

(n) 

suppose that Xq ^ x as n — )■ 00. Then, up to a time-rescaling depending on 
A and the scaling of the stable process above, the processes X^^^ converge 
to a Markov process X that is a critical (1 + /3)-stable continuous state 
branching process, and whose distribution is determined by the Laplace 
transform 

(1.1) E[e-''^*jXo = 2;] = expf 



(l+0/3t)l//3 

If /3 = 1, this is Feller's critical continuous state branching process [17, 22]. 
(Note that time here is scaled by a factor of 2 relative to some other authors 
so that the generator of our "Feller continuous state branching process" is 

Let T = inf {t > : = 0} denote the extinction time of X (it is not hard 
to show that = for all t>T). Taking ^ ^ 00 in (1.1) gives 



[t > t\Xo = x} = 1 — exp 



X 

"t^ 



so X dies out almost surely. However, it is possible to condition X to live 
forever in the following sense: 

hm E[f{Xt)\Xo = x,T>T] = -E[f{Xt)Xt\Xo = x]. 

T^oo X 

Thus if Pt{x' ,dx") are the transition probabilities of X, then there is a 
Markov process Y with transition probabilities 

Qtiy',dy") = ^Ptiy',y")y". 

y 

The process Y is the critical (1 + /3)-stable continuous state branching pro- 
cess X conditioned on nonextinction. The distribution of Yt is determined 
by its Laplace transform 

(1.2) n^M-OYt)\Y, = y] = exp(^-^^^^^^^) (1 + 16^)-^^+')/^ 

(see Section 6). Moreover, it is possible to start the process Y from the 
initial state Iq = 0, and the formula (1.2) continues to hold for y = 0. The 
super-process generalization of this construction was considered for /3 = 1 in 
[11, 14, 15] and for general /3 in [16]. 

For /3 = 1, the conditioned process Y can be described informally as a sin- 
gle "immortal particle" constantly throwing off infinitesimally small masses 
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with each mass then evolving according to the dynamics of the uncondi- 
tioned process. These infinitesimal masses can be interpreted as the single 
progenitors of families whose lineage splits from the immortal particle at 
the birth time of the progenitor and are eventually doomed to extinction. 
Most such families die immediately, but a rare few live for a noninfinitesimal 
amount of time. More formally, there is a cx-finite measure v on the space of 
continuous positive excursion paths 

5° := {n G C(M+, M+) : uq = & 37 > s.t. tit > ^ < t < 7}, 

such that if n is a Poisson point process on IR+ x with intensity \ ®v 
where A is Lebesgue measure, and (Xj)t>o is an independent copy of X 
begun at = y, then the process 




has the same distribution as (Yi)f>o begun at Yq = y (see [15]). A point 
(s,ti) S n corresponds to a family that grows to nonnegligible size; the time 
s records the moment the family splits off from the immortal particle, and the 
value Ur of the trajectory u gives the size of the family r units of time after it 
split off. The family becomes extinct after the period of time ^{u) := inf{r > 
0:it(t) = 0,Vt > r}. The a-finite measure v may be identified explicitly, but 
it suffices to remark here that it is Markovian with transition probabilities 
the same as those of the unconditioned process X — in other words, v arises 
from a family of entrance laws for the semigroup of X. The process (Xf)4>o 
records the population mass due to descendants of individuals other than 
the immortal particle who are present at time 0. 

An analogous description of the conditioned process Y for the case /3 € 
(0,1) is presented in [16]. There is again a single immortal lineage, but 
now families split off from that lineage with a noninfinitesimal initial size, 
refiecting the heavy-tailed offspring distributions underlying these models. 
More precisely, a decomposition similar to (1.3) holds, but the Poisson point 
process 11 is now on M4. x E where E = {u^ Z)(M+, M_|_) : 37 > s.t. > 44> 
< t < 7}, the set of cadlag paths starting above zero that eventually hit 
zero. The nondecreasing process (Mj)t>o where 

Mt := ^ uo 

{s,u)enn[o,i]x£: 

is the total of the initial family sizes that split off from the immortal particle 
in the time interval [0,t]. It is a stable subordinator of index /3. 

Suppose now that /3 € (0, 1] is arbitrary. Take Iq = so that = in 
the decomposition (1.3) and all "individuals" belong to families that split 
off from the immortal particle at times s > 0. Extend the definition of 7(n) 
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At 



Fig. 1. The points of A are marked "x"; the sample path of the process A is drawn with 
a solid line; the oldest extant family at time t is represented by the point in A labeled "O" 
and the left-leaning wedge A{t,0) is the darkly shaded region with apex at (t,0). At time 
T, the family represented by O will die out and the family represented by N will contain 
the MRCA because O is on the boundary of (the lightly shaded) A(T, 0). If the coordinates 
of O are {u,y) and the coordinates of N are {v,z), then the size of the jump that A makes 
at T is At- At- = u - v = {T - v) ~ [T - u) = [T ~ v) ~ y . 

given above for u € £"0 to u € in the obvious way. The individuals, besides 
the immortal particle alive at time t > 0, belong to families that correspond 
to the subset At ■= {{s,u) Gn:0<t — s< jiu)} of the random set 11. At 
time t, the amount of time since the most recent common ancestor of the 
entire population lived is At := sup{i — s : (s, n) € At}- As depicted in Figure 
1, the MRCA age process {At)t>o has saw-tooth sample paths that drift up 
with slope 1 until the current oldest family is extinguished, at which time 
they jump downward to the age of the next-oldest family. 

It is not necessary to know the Poisson point process 11 in order to con- 
struct the MRCA age process {At)t>o- Clearly, it is enough to know the 
point process A on M+ x M+4- given by A := {{s,'y{u)) : S 11}. Indeed, 
if we define the left-leaning wedge with apex at {t,x) to be the set 

(1.4) A(t,x) ■.= {{u,v) eR^:u<t &i u + v>t + x}, 

then 

At = t- inf {s : 3x > s.t. (s, x) G A n A{t, 0)} 

(see Figure 1). 

Note that A is a Poisson point process with intensity A ® /i where fi is the 
push-forward of u by 7; that is, //((i,oo)) = ;/({n:7(n) > t}). We will show 
in Section 5 that /i((t,oo)) = (1 + /3)/(/?t). 

With these observations in mind, we see that if A is now any Poisson point 
process on M-|_ x M-i-i- with intensity A (8) ;U where fi is any measure on M_|__|_ 
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with /i(M4._|„) = oo and < oo)) < oo for all x > 0, then the construction 
that built {At)t>o from the particular point process A considered above will 
still apply and produce an M+-valued process with saw-tooth sample paths. 
We are therefore led to the following general definition. 

Definition 1.1. Let A be a Poisson point process on M+ x M++ with 
intensity measure X0 fi where A is Lebesgue measure and is a c-finite 
measure on with ^(M_(_-|_) = oo and oo)) < oo for all x > 0. Define 

{At)teR+ by 

At:=t- inf{s > : 3x > s.t. (s, x) G A n A(t, 0)}, 
where A(t,0) is defined by (1.4), and = if A n A(t,0) is empty. 

We will suppose from now on that we are in this general situation unless 
we specify otherwise. We will continue to use terminology appropriate for 
the genealogical setting and refer to {At)t>o as the MRCA age process and fi 
as the lifetime measure. We will assume for convenience that the measure fj, 
is absolutely continuous with a density m with respect to Lebesgue measure 
that is positive Lebesgue almost everywhere. It is straightforward to remove 
these assumptions. 

The strong Markov property of the Poisson point processes A implies that 
{At)t>o is a time- homogeneous strong Markov process. In particular, there 
is a family of probability distributions (P^)a:eK_^ on the space of M+-valued 
cadlag paths with interpreted in the usual way as the "distribution of 
{At)t>o started from Aq = x." More concretely, the probability measure is 
the distribution of the process {Af)t>o defined as follows. Let A^ be a point 
process on [— x, oo) x M+4- that has the distribution of the random point 
set {{t — x,y): {t,y) G A} U {(— x,Z)} where Z is an independent random 
variable that is defined on the same probability space {Q,J^,¥) as the point 
process A, takes values in the interval (x, 00) and has distribution 

F{Z<z} = ^i{{x,z])/fi{{x,^)). 

Then we can define 

Af :=t- inf {s > -x : 3y > s.t. (s, y) G A"^ n A(t, 0)}, t > 0, 

where we adopt the convention that A^ : = if the set above is empty [see the 
proof of part (a) of Theorem 1.1 below for further information]. From now 
on, when we speak of the process {At)t>o we will be referring either to the 
process constructed as in Definition 1.1 from the Poisson process A defined 
on some abstract probability space {Q,J-,¥) or to the canonical process on 
the space of cadlag M+-valued paths equipped with the family of probability 
measures (P^)a;>o- This should cause no confusion. 

We note in passing that the analogue of {At)t>o in the constant population 
size Wright-Fisher setting is not Markov (see Remark 4.1.3 of [26]). 

We prove the following properties of the process {At)t>o in Section 2. 
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Theorem 1.1. (a) The transition probabilities of the time-homogeneous 
Markov process {At)t>o have an absolutely continuous part 

F-{At G dy} = ^^^-^'"^^.y exp (- f^' ^((z, oo)) dz\ /.((y, c^)) dy, 

for y < X + t and a single atom 

F^A = x + t}= ^%-^'^'^^\ 
/i((x,oo)) 

(b) The total rate at which the process {At)t>o jumps from state x > is 

m{x) 
n{{x,oo)y 

and when the process jumps from state x > 0, the distribution of the state 
to which it jumps is absolutely continuous with density 

y^exp(^-j n{{z,oo)) dz^fi{{y, oo)), 0<y<x. 

(c) The probability F^{3t > : = 0} that the process {At)t>o returns to 
the state zero is positive if and only if 

J exp^y ix{{y, oo)) dy^ dx < oo. 

(d) // 

expf-y n{{y,oo))dy] dx = oo, 

then for each x > the set {t >0:At = x} is ¥^ -almost surely un- 
bounded. Otherwise, limt-^c>o At = oo, ¥^ -almost surely, for all x > 0. 

(e) A stationary distribution vr exists for the process (At)t>o if o-'^d only if 



oo 



fj,{{z, oo)) dz < oo, 
in which case it is unique, and 

Tr{dx) = /u((rE, oo)) exp ^— J fJ^dz, oo)) dz ) dx. 
(f) If {At)t>o has a stationary distribution n, then 

dTv(lP"{^t G •},vr) < 1 -expj^-^^ ^l{iy , oo)) dy^ 



n{{y,oo))dy x 



fx{[x,oo)) 

where d^y denotes the total variation distance. In particular, the distri- 
bution of At under converges to vr in total variation as t—^oo. 
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Specializing Theorem 1.1 to the case when A is the MRCA age process 
of the conditioned critical (1 + /3)-stable continuous state branching pro- 
cess gives parts (a) to (d) of the following result. Part (e) follows from an 
observation that a space-time rescaling of this MRCA age process is a time- 
homogeneous Markov process that arises from another Poisson process by 
the general MRCA age construction of Definition 1.1. The proof is in Section 
5. 



Corollary 1.1. Suppose that A is the MRCA age process associated 
with the critical (1 + (3) -stable continuous state branching process. 

(a) The transition probabilities of the process A have an absolutely contin- 
uous part 

and a single atom 

F^At = x + t}- 



X -\-t 

(b) The total rate at which the process A jumps from the state x > is 1/x, 
and when it jumps from state x> 0, the distribution of the state to which 
it jumps is absolutely continuous with density 

+ 0<y<x. 

(c) The probability W^{3t > : = 0} that the process A returns to the state 
zero is 0. 

(d) For each x>0, liuit^oo At = oo, Y'^ -almost surely. 

(e) The process 

indexed by the whole real line is a time-homogeneous Markov process 
under for any x>0, and it is stationary when x = 0. Moreover, At/t 
converges in distribution to the Beta(l -|- 1//3, 1) distribution as t oo 
under P^' for any x >0, and At/t has the Beta(l -|- 1//5, 1) distribution 
for all t> when x = 0. 

Note that the sample paths of {At)t>o have local "peaks" immediately 
before jumps and local "troughs" immediately after. We investigate the dis- 
crete time Markov chain of successive pairs of peaks and troughs in Section 
4. We also consider the jump heights and inter-jump intervals and describe 
an interesting duality between these sequences in Section 3. 
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Finally, recall that the Bessel-squared process in dimension 7, where 7 
is an arbitrary nonnegative real number, is the M+-valued diffusion process 
with infinitesimal generator 1x cP' /dx^ + '^d/dx. When 7 is a positive integer, 
such a process has the same distribution as the square of the Euclidean norm 
of a Brownian motion in R'^. Feller's critical continuous state branching 
process is thus the zero-dimensional Bessel-squared process, modulo a choice 
of scale in time or space. It was shown in Example 3.5 of [28] that for < 7 < 
2, the Bessel-squared process with dimension 7 conditioned on never hitting 
zero is the Bessel-squared process with dimension 4 — 7. Thus for /3 = 1, the 
conditioned process Y is the four-dimensional Bessel-squared process. We 
introduce a new family of processes in Section 6 that are also indexed by a 
nonnegative real parameter and play the role of the Bessel-squared family 
for values of /? other than 1. These processes will be studied further in a 
forthcoming paper. 

We end this introduction by commenting on the connections with previous 
work. First, we may think of each point (s,x) € A as a "job" that enters a 
queue with infinitely many servers at time s and requires an amount of time 
X to complete. We thus have a classical M/G /oo queue [31], except we are 
assuming that the total arrival rate of jobs of all kinds is infinite. With this 
interpretation, the quantity At is how long the oldest job at time t has been 
in the queue. Properties of such M/G/ 00 queues with infinite arrival rates 
have been studied (see, for example, [10]), but the age of the oldest job does 
not appear to have been studied in this context. 

Second, note that the process {At)t>o is an example of a piecewise de- 
terministic Markov process: it consists of deterministic flows punctuated by 
random jumps. Such processes were introduced in [5] and studied further 
in [6] (see also [20], where the nomenclature jumping Markov processes is 
used). The general properties of such processes have been studied further 
in, for example, [3, 4, 7]. 

Last, piecewise deterministic Markov processes like {At)t>o that have pe- 
riods of linear increase interspersed with random jumps have been used to 
model many phenomena, such as stress in an earthquake zone [2], conges- 
tion in a data transmission network [8] and growth-collapse [1]. They also 
have appeared in the study of the additive coalescent [12] and M-tree- valued 
Markov processes [13]. 

2. Proof of Theorem 1.1. (a) Suppose that {At)t>o is constructed from 
A as in Definition 1.1. For s>x consider the conditional distribution of Ag-^-t 
given As = X. The condition As = x is equivalent to the requirements that 
there is a point (s — x, Z) in A for some Z > x and that, furthermore, there 
are no points of A in the left-leaning wedge A(s — x,x). The conditional 
probability of the event {Z > z}, given As = x, is fi{{z, 00)) / fi{{x, 00)) for 
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z > X. If Z > X + t, then Ag+t = x + t. Otherwise, Ag^t < x + t. The second 
claim of part (a) follows immediately. 

Now consider F{As^t € dy\As = x} ior y < x + t. This case is depicted in 
Figure 2. By construction, A^^t = 2/ if and only if there is a point {s + t — 
y, W) £ A for some W > y, and there are no points of A in A(s + t — y,y). 
From above, the condition Ag = x requires there to be no points of A in 
the wedge A(s — x,x) (the lightly shaded region in Figure 2). Therefore, if 
As = X, then As+t = ?/ if and only if Z <x + t, there are no points of A in 
the darkly shaded region of Figure 2 and there is a point of A of the form 
{s + t — y,w) with w > y. 

Now, the conditional probability of the even {Z <x + t} given As = x is 



The probability that no points of A are in the darkly shaded region is 



The probability that A has a point in the infinitesimal region [s + t — y,s + 
t — y + dy] X {y,oo) is iJ,{{y , oo)) dy . Multiplying these three probabilities 
together gives the first claim of part (a). 

(b) Both claims follow readily by differentiating the formulae in (a) at 
t = 0. It is also possible to argue directly from the representation in terms 
of A. 



H{{x,x + t]) 

//((X,CX))) 




given no point point must fall 

in this region in this interval 



no point can fall 
in this region 




point must fall 
in this half hne 



X 



s 



y 



s+t 



Fig. 2. The computation off{As+t G dy\As — x}. 
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(c) Suppose that {At)t>o is constructed from A as in Definition 1.1. Note 
that 

{t>0:^t = 0} = R++\ IJ {t,t + x). 

The question of when such a "Poisson cut out" random set is almost surely 
empty was asked in [23] , and the necessary and sufficient condition presented 
in part (c) is simply the one found in [30] (see also [18]). 

(d) For X > 0, the random set {t >0:At = x} is a discrete regenerative 
set under that is not almost surely equal to {0}. Hence this set is a 
renewal process with an inter-arrival time distribution that possibly places 
some mass at infinity; in which case, the number of arrivals is almost surely 
finite with a geometric distribution, and the set is almost surely bounded. 

Suppose that the set {t >0:At = x} is P^-almost surely unbounded for 
some X > 0. Let = Tq < Ti < • • • be the successive visits to x. It is clear that 
P^{3t € [0, Ti] ■.At = y}>0 for any choice of y > and hence, by the strong 
Markov property, the set {t >0 : At = y} is also P^-almost surely unbounded. 
Another application of the strong Markov property establishes that the set 
{t >0: At = y} is P^-almost surely unbounded. Thus the set {t>0:At = x} 
is either unbounded P^-almost surely for all x > or bounded P^-almost 
surely for all x > 0. 

Recall that away from the set {t>0:At = 0} the sample paths of A are 
piecewise linear with slope 1. It follows from the coarea formula (see, e.g., 
Section 3.8 of [24]) that 

roo /"OO 

/ f{At)dt= f{y)i^{t>Q:At = y]dy 
Jo Jo 

for a Borel function / : IR+ — )• M-j.. Hence by Fubini's theorem, 

fivWmt >0:At = y}]dy = f{y) ^"^^^^^'^V^ dt dy 

for any x > 0. It follows from the continuity of the transition probability 
densities that 



W[#{t>Q:At = y}]= f 

Jo 



^'{Atedy}^^ 



dy 



for all x,y > and, in particular, that the expected number of returns to 
X > under P^ is 







dx 



Using the expression from part (a) and the argument above, this quantity 
is infinite, and hence the number of visits is P^-almost surely infinite, if and 
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only if 

exp( — / fi{{y, oo)) dy ] du = oo. 



If the set {t>0:At = x} is P^-almost surely bounded for all x > 0, then 
by an argument similar to the above, the set {t >0:At = y} is P^-almost 
surely bounded for all x,y > 0. It follows that, for all x > 0, P^-almost surely 
all of the sets {t>0:At = y} are finite. This implies that liuit^oo At exists 
P^-almost surely, and the limit takes values in the set {0,oo}. However, it is 
clear from the Poisson process construction that does not occur as a limit 
with positive probability. 

(e) Suppose there exists a probability measure vr on such that 

/ F{Atedy\Ao = x}7ridx)=7r{dy), yGM+. 
Jr+ 

Taking t — t- oo in part (a) gives 

7r(dy)=lim f ^«-'- + *]) 



rx+t 

xexpl— / fj.{{u, oo)) du ]Tr{dx)fi{{y,oo)) dy 



y 



/•oo 

0, if / ^{{u, oo)) du = oo, 
Jy 



exp ( - y oo)) du ] fi{{y, oo)) dy, 

otherwise. 



Therefore, a stationary probability distribution exists if and only if J^°°/i(( 



u. 



oo)) du < oo, and if a stationary distribution exists, then it is unique. 

(f) It will be useful to begin with a concrete construction of a stationary 
version of the process A in terms of a Poisson point process. Suppose that 
fi{{u,oo)) du < oo for all x > 0, so that a stationary distribution exists. 
Let be a Poisson point process on M x M++ with intensity measure A (8)^. 
Define {Af)tm by 

:= t - inf{s : 3x > s.t. (s, x) G n A{t, 0)}. 

The condition on /i ensures that almost surely any wedge A(t,x) with x > 
will contain only finitely many points of A'*"'', and so {A'j~^)t£R is well defined. 
The process (A^)jg]R is stationary and Markovian, with the same transition 
probabilities as {At)t>o- 

Recall the construction of the process A^ started at x for x > that was 
described preceding the statement of Theorem 1.1. Construct the point pro- 
cess A^ that appears there by setting A^ := {{t, y) G A^ : t > — x}U{(— x, Z)} 
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where Z is an independent random variable with values in the interval (x, oo) 
and distribution P{Z > z} = ^((z, oo))/ ii{{x, oo)). 

By construction, Af = A^^ for all i > T where T is the death time of all 
families alive at time —x in either process: 

T := mf{t > : Z < t + X and n A(-x, t + x) = 0}. 

Thus 

dTy(F^{Ate-},7r)<¥{Af^Ar} 
< F{T > t} 

= 1 - F{Z < t + x}P{A^ n A(-x, t + x) = 0} 

and part (f) follows. 

3. Duality and time-reversal. Suppose in this section that J°° fj,{{y, 
oo)) dy < oo for all x > 0, so that, by part (e) of Theorem 1.1, the pro- 
cess A has a stationary distribution. Let {Af^)t^M. be the stationary Markov 
process with the transition probabilities of A that was constructed from the 
Poisson point process in the proof of part (f) of Theorem 1.1. 

Define the dual process (A^)j6ir by A^ := inf{s > : A{t, s) n A^ = 0}. 
See Figures 3 and 4. Thus A'^ is the amount of time that must elapse after 
time t until all families alive at time t have died out, or, equivalently, until 
the MRCA for the population lived at some time after t. The cadlag Re- 
valued process (^^)jgiR has saw-tooth sample paths that drift down with 
slope —1 between upward jumps. 

Proposition 3.1. The dual process, (A^)tgR, has the same distribution 
as the time-reversed process, (A^)tgK, where At:=liniuit^t^u- 

Proof. Define bijections ip and o" of M x M++ by ip{t,x) := {t + x,x) 
and a{t,x) := (— t,x). In the usual manner, we may also think of ip and a 
as mapping subsets of M x M++ to other subsets of M x M+_|_. Note that ip 
maps left-leaning wedges to right-leaning wedges, and a maps right-leaning 
wedges to left-leaning wedges. Thus the composition a oip maps left-leaning 
wedges to left-leaning wedges. More precisely, 

a o {p[A{t, x)) 

= {a oip{s,y):{s,y) G A(t,x)} 

(3.1) ={{-{s + y),y):s<t &i s + y>t + x} 

= {{u,v):u <-{t + x) k u + V > -t = -{t + x) + x} 

= A{aoip{t,x)), 
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Fig. 3. The process A'^ and the dual process . The points in the shaded area represent 
the families alive at time t, and the solid line is the sample path of At* . The point marked 
"O" is the oldest living family at time t; the point marked "Y" is the family extant at time 
t that will live the longest into the future. 




Fig. 4. The (coupled) processes A"^ and A*^ for the same set of family lifetimes (the 
"x"s). The paths of A^ begin at points in r(A") and descend; while the paths of A*^ 
ascend to points m (fi{r{A^)) . The mapping ip is shown by the horizontal dotted lines in 
the lower diagram. 
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where we stress that A{a o ip[t,x)) is the left-leaning wedge with its apex 
at the point a o (f[t,x) € M x M4.4.. Define a map T that takes a subset of 
R X M_|_+ and returns another such subset by 

r{s) := {{t,x) e S : A{t,x) n s = 0}. 

The points of r(A^) correspond precisely to those families that at some 
time will be the oldest surviving family in the population. These points 
determine the jumps of both the MRCA process and the dual process: as 
can be seen with the help of Figure 4, the linear segments of the paths of 
the dual process {Al*)t£M each begin at a point in r(A^) and descend with 
slope —1 whereas the linear segments of the paths of the MRCA process 
{AP')t£M. ascend with slope +1 to points in if oT{A^). 

This implies that the path of the time-reversed process, {At)t£M., begin 
at points in a o o r(A^) and descend with slope —1, and so the points 
of 0" o o r(A^) determine the path of {At)t£R in the same way that the 
points of r(A^) determine the path of {At)t£K- Therefore, all we need to do 
is show that a o (po r(A^) has the same distribution as r(A'^). 

To see this, first note that, for an arbitrary subset 5* C M x M_|_+, 

a o Lpo r(5) 

= {(To ip{t, x) : {t, x)eS and A(t, x) H S = 0} 
= {{t,x):ip~^ oa~^{t,x) eS and A{ip~'^ o a'"'^{t,x)) D S = 0} 
= {{t,x):ip'~^ oa~^(t,x) e S and 99"^ o cr~^(A(t, x)) n S = 0} 
= {{t, x) :{t,x) £ao ip{S) and A{t,x)nao ip{S) = 0} 

= r o (7 o ip{s), 

where the third equality follows from identity (3.1). Thus a o ip o r(A^) = 
r o CT o {p(A'^). However, both maps and a preserve the measure A (8) ^ and 
hence a o ip[A^) has the same distribution as A'*^. Therefore, a oipo r(A^) 
has the same distribution as r(A^), as required. □ 

Remark 3.1. There is an interesting connection between the jump sizes 
and the inter-jump intervals, stemming from the observation that the paths 
of A^ and A^ have the same sequences of "trough" and "peak" heights, 
while the roles of the jump sizes and inter-jump intervals for the two are 
exchanged. The following situation is depicted in Figure 4. To explain the 
connection, suppose that T G {t € M:^^ 7^ ^t^} is a jump time for the 
process A^. Let T< := sup{t <T:Alt^ A^^} and T> := mf{t >T:A^j^ 
A^*} be the jump times on either side of T. Put L := Ajt_ and R := A^ , 
and define L^, R^, and B?' as the analogous left limits and values of 
A*^ at the times and T^. Write A := L — R for the size of the jump 
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at time T and H := T — for the length of the time interval since the 
previous jump. Observe that T — L is a jump time for the dual process A'^ 
with A'^_j^^_ = i?< and A'^_j^^ = L. Moreover, 

H = L - = - 

and 

A = (r> - L>) -{T-L) = mf{t >T-L:Atl^ If} - {T - L). 

Note that the map T ^ T — L sets up a monotone bijection between the 
jump times of the process A^ and those of the process A'^ . It thus follows 
from Proposition 3.1 that the point processes, 

{(T, r - r<, - A^) : A^_ + 

and 

{(T, A^_-A1^,T- T<) : At;_ ^ A^}, 
have the same distribution. 

4. Jump chains. Suppose again that f^fi{{y,oo))dy < oo for all x > 
so that the process A has a stationary distribution. Recall the stationary 
Markov process, {At*)t£R, with the transition probabilities of A that was 
constructed from the Poisson point process, A^, in the proof of part (f) of 
Theorem 1.1. 

For t G M, denote by Jt := inf{u > 0:^1^ ^ ^tt-} the next jump time 
of A'^ after time t. Define an increasing sequence of random times, < 
Tq <Ti < ■ ■ ■, by To := Jo and Tn+i := Jt„ for n > 0. Put L„ := A^^^_ and 
Rn:=A^^. Thus the sequences {Ln)^=Q and (i?n)^o record, respectively, 
the "peaks" and the "troughs" of the path of A"^ that occur between the 
times and sup„T„. 

The next result can be proved along the same lines as part (a) of Theorem 
1.1, and we leave the proof to the reader. 

Proposition 4.1. The sequence, {Lq, Rq, Li, Ri, . . .), is Markovian with 
the following transition probabilities: 

¥{Rnedy\Ln = x} = iJ,{{y,oo))exp^- j ij,{{u, oo)) du^ dy, 0<y<x, 



and 



TThiz) 

^{Ln+i^dz\Rn=y} = — —dz, z>y. 



In particular, the sequence of pairs i?„))^Q is a time-homogeneous 

Markov chain. 
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Now we may compute the transition probabilities of the peak and trough 
chains. By Proposition 4.1, 

¥{Ln+iedz\Ln = x}/dz 

/ \ m(z)l < 

(4.1) =J^ exp^-y^ f,{iu,oc))dujf,{{y,oc))^^^^dy 

f-xAz / i-x \ 

= m{z) / expl — / fx{{u, oo)) duj dy 



and 

¥{Rn+iedz\Rn = x}/dz 



y 



(4.2) =/ — yy/x((2;,oo))exp( - / fi{{u,oo))du]ly>^dy 



m{y) 



y 



= n{{z,oo)) —- --exp - / fj,{{u,oo))du] dy. 

It follows from (4.1) that the peak chain, (in)^O' X-irreducible where 
A is Lebesgue measure on That is, if A is a Borel subset of with 

A(^) > 0, then, for any x € ^++, there is positive probability that the the 
peak chain begun at x will hit A at some positive time (see Chapter 4 of 
[25] for more about this notion of irreducibility) . It follows that the peak 
chain is either recurrent, in the sense that 



J2^{Ln^A\Lo = x} = oo 



n=0 



for all X € M++ and all Borel subsets of A C M++ with A(^) > 0, or it is tran- 
sient, in the sense that there is a countable collection of Borel sets {Aj)°^-^ 
and finite constants {Mj)^^-^ such that Uj^i^j =I^++ ™d 



x& 



00 

sup V F{Ln G Aj\Lo = x}< Mj 

*++ n=0 

(see Theorem 8.0.1 of [25]). 

The peak chain is strong Feller; that is, the function, 

xi->E[/(L„+i)|L„ = a;], 

is continuous for any bounded Borel function /. Also, because the support 
of A is all of M_|_+, if the peak chain is recurrent, then each point x of M_|__(_ 
is topologically recurrent in the sense that 



^P{L„ G U\Lo = x} = oo 



n=0 
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for every open neighborhood U of x. Hence, by Theorem 9.3.6 of [25], if the 
peak chain is recurrent, then it is Harris recurrent, which means that given 
any Borel set A with X{A) > 0, the chain visits A infinitely often almost 
surely starting from any x. Moreover, the chain is recurrent (equivalently, 
Harris recurrent) if and only if it is nonevanescent; that is, started from any 
X, there is zero probability that the chain converges to or oo (see Theorem 
9.2.2 of [25]). 

If the peak chain is recurrent (equivalently, Harris recurrent or nonevanes- 
cent), then it has an invariant measure that is unique up to constant mul- 
tiples (see Theorem 10.4.4 of [25]). If the invariant measure has finite mass, 
so that it can be normalized to be a probability measure, then the chain is 
said to be positive, otherwise the chain is said to be null. 

Conversely, if the peak chain has an invariant probability measure, then 
it is recurrent (equivalently, Harris recurrent or nonevanescent) (see Propo- 
sition 10.1.1 of [25]). 

All of the remarks we have just made for the peak chain apply equally to 
the trough chain (i?„)J^g. Recall that we are in the situation when A has a 
stationary version, so the transience or recurrence of L and R depends on 
their behavior near zero. 

Proposition 4.2. Consider the two Markov chains, (L„)^q and (i?ri)5JLo- 

(a) Both chains are transient if and only if 

exp^y fi{{y, oo)) dy^ dx < (X). 

(b) Both chains are positive recurrent if and only if 
m(x)exp(— / fi{{y, oo)) dy \ dx < oo. 







(c) Both chains are null recurrent if and only if both 
and 

m(x)exp(— / fi({y, oo)) dy \ dx = oo. 



I 

Jo 



Proof. Consider the set Z := {teR: A^ = 0}. It foUows from part (c) 
of Theorem 1.1 that ¥{Z ^ 0} > if and only if 



(4.3) 



y exp^y fi{y , oo) dy^ dx < oo . 
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By the stationarity of {Al*)t^^ and the nature of its sample paths, it is 
clear that for x > 0, the set {t € R : A^* = x} is unbounded above and below 
almost surely [this also follows from part (d) of Theorem 1.1]. It follows from 
a simple renewal argument that if (4.3) holds, then Z is unbounded above 
and below almost surely. 

Because the paths of (A^)tgK increase with slope 1 in the intervals 
[T„,T„+i), it follows that if condition (4.3) holds, then lim„__s.oo = inf{t > 
: A'^ = 0} < oo almost surely, and lim„^oo Ln = lim,„^oo Rn = almost 
surely. In this case, both chains are evanescent, and hence transient. 

On the other hand, if (4.3) does not hold, then lim„_j,ooTk = oo. More- 
over, the set {t € M : A^^ = x} is almost surely unbounded above and below 
for any x > 0, as we observed above. If we split the path of {A'^)t,^K into ex- 
cursions away from x, then each excursion interval will contain only finitely 
many jumps almost surely and, because the excursions are independent and 
identically distributed, it cannot be the case that L„ or Rn converges to 
or oo with positive probability. Thus, both chains are nonevanescent and 
hence recurrent. 

It is clear from (4.1) that the kernel giving the transition densities of 
the peak chain (-Ln)5JLo self-adjoint with respect to the measure having 
density 

p{x) = m{x) ex.pl — / /i((ii, oo)) dn 



with respect to Lebesgue measure, and so this measure is invariant for the 
peak chain. Clearly, p{x) dx < oo if and only if the condition in part (b) 
holds in which case the peak chain is positive recurrent. Otherwise, the peak 
chain is either null recurrent or transient, and so part (a) shows that the 
peak chain is null recurrent if the two conditions in part (c) hold. 

Similarly, it is clear from (4.2) that the kernel giving the transition densi- 
ties of the trough chain (i?„)^Q is self-adjoint with respect to the measure 
having density 

g(x) = //((x, oo))^ exp^— y fj,{{u,oo))du 

with respect to Lebesgue measure, and so this measure is invariant for the 
trough chain. An integration by parts shows that q{x) dx < oo if and only 
if the condition in part (b) holds, and so the trough chain is positive if and 
only if the peak chain is positive. Alternatively, we can simply observe from 
Proposition 4.1 that integrating the conditional probability kernel of Rn 
given Ln against an invariant probability measure for the peak chain gives 
an invariant measure for the trough chain, and integrating the conditional 
probability kernel of Ln+i given Rn against an invariant probability measure 
for the trough chain chain gives an invariant measure for the trough chain, 
so that one chain is positive recurrent if and only if the other is. □ 
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Remark 4.1. If m{x) = ax ^ for x S (0,1], then both the peak and 
trough chains are: 

(1) transient <^ < a < 1; 

(2) nuh recurrent <4> a = 1; 

(3) positive recurrent <^ a > 1. 



Remark 4.2. It follows from parts (b) and (e) of Theorem 1.1 that the 
stationary point process {t € M : A^^ ^ has intensity 

p ■= m(a;)exp(— / p,{{u, oo)) du \ dx 



and so the peak and trough chains are positive recurrent if and only if p is 
finite. Suppose that p is finite and consider the point process 

E := {{t, Atl,At') G R X R+ X M+ : / }. 

The companion Palm point process T has its distribution defined by 

'-{n:0<r„<l} 

where O^B = {{t - s, i, r) : {t, i, r) G 5} for 5 C M x M+ x M+. Enumerate the 
points of T as ((r„, L„, Rn))nez where • • • < T_i < Tq = < Ti < • • •. A fun- 
damental result of Palm theory for stationary point processes says that the 
random sequence ((T„ — T„_i, L„, is stationary and that the distri- 

bution of the point process H may be reconstructed from the distribution of 
this sequence (see, e.g.. Theorem 12. 3. II of [9] or [21]). It is clear that the 
stationary random sequences (Lri)5^o (-Rn)^o have the same distribu- 
tion as the peak and trough chains started in their respective stationary 
distributions. 



5. The (1 -|- /3)-stable MRCA process. In this section we specialize to 
the motivating example of the MRCA process of a critical (1 -1- /3)-stable con- 
tinuous state branching process conditioned to live forever. Recall that the 
unconditioned continuous state branching process has Laplace transforms 
(1.1), and the conditioned process has Laplace transforms (1.2). For /3 = 1, 
the unconditioned process has generator x with this choice of time scale. 

As mentioned in the Introduction, the set of points {t^x) G M_|_ x 
where t is the time that a family splits from the immortal lineage and x is 
its total lifetime, is a Poisson point process with intensity measure \ ® jjb for 
some (T-finite measure p. 
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Lemma 5.1. The lifetime measure fj, associated with the critical (1 + /3)- 
branching process conditioned on nonextinction is given by 

1 + /3 

/u((x,oo)) = x>0. 

Proof. As we remarked in the Introduction, 
(5.1) F{Xt>0\Xo = x} = l-exp\ 



t^lP 



First consider the case of /? = 1. Recall from the Introduction that if 11 is 
a Poisson point process on R_|„ x £^ with intensity A (8) z^, then 

(5-2) ( J] U(t_^)vo 

has the same distribution as the conditioned process, (l^)t>o, with Yq = 0, 
and recall that ^ is the push- forward of v by the total lifetime function 7. 
Also, 

(5.3) ( E ^* 

has the same distribution as the unconditioned process [Xt)t>Q with Xq = y 
(see [15]). The factor of 2 differs from [15] and arises from our choice of time 
scale. Therefore, 

W{Xt > 0|Xo = y]= P{3(s, li) G n : s < y/2 and 7(u) > t} 

= 1 -exp(-y^((t,oo))/2), 

and comparing with (5.1) gives //((t,oo)) = 2/t. 

Now take /3 G (0,1). It is shown in [16] that the mass thrown off the 
immortal lineage is determined by the jumps of a stable subordinator: if Mg 
is the amount of mass thrown off during the time interval [0, s] , then 



E[e~^*'^^] =exp( -s 



1 + 



/5 



where v{dx) = y\i-i3) dx is the Levy measure of the subordinator. 

Since the jump size of the subordinator corresponds to the initial size of 
the new family, the lifetime measure /i is given by 



^l{{t,oo)) 



poo 

/ F{Xt >0\Xo = x} i^{dx) 
Jo 
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and so, from the above and an integration by parts, 

fi{{t, oo))= / (l-e-^/*'^'')z^((ix) 



1/^, 

(i-e 

^0 

_ l + f3 

~ pt ■ □ 

Proof of Corollary 1.1. Parts (a) to (d) follow immediately from 
Theorem 1.1. Part (e) will also follow from parts (e) and (f) of Theorem 1.1 
after the following time and space change. 

Define a new time parameter n by t = e". If the MRCA at time t lived 
at time t — x on the original scale, then on the new time scale she lived at 
time u — y where t — x = e"~^. Solving for y, the MRCA age process in the 
new time scale is the process {Bt)t>o given by Bu = — log(l — e~"^e")- The 
process {Bt)t>o is obtained by applying the construction (1.1) to the point 
process given by 

logs,logfl + ^^^,(s,2;) G a|, 

which is a Poisson point process on M x M_|_ with intensity measure A p 
where 

1 + /3 

p{{y,oo)) = ^^^y y£R++. 

Note that, in general, a time and space change of the Poisson process A 
gives a new Poisson point process, but the resulting intensity measure will 
not typically be of the form A (8) k for some measure k: it is a special feature of 
fi and the transformation that the product measure structure is maintained 
in this case. 

It is straightforward to check parts (e) and (f) of Theorem 1.1 that {Bt)t>o 
has the stationary distribution 

Tr{dx) = i±^e-^(l - e"^)^/^ dx 

and that the distribution of Bf converges to vr in total variation as t — ?> oo. 

Part (e) of the corollary then follows from the observation that ^ = 1 — 
-a 



e '°s(t) and an elementary change of variables. □ 



6. An analogue of the Bessel-squared family. In this last section, we de- 
termine the Laplace transform of the (1 + /3)-stable continuous state branch- 
ing process conditioned on nonextinction. For /3 = 1, the unconditioned pro- 
cess is the Bessel-squared process with dimension and it is well known 
that the conditioned process is the Bessel-squared process with dimension 
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4. A comparison of Laplace transforms will suggest that the unconditioned 
and conditioned processes for a given general (3 may also be embedded in a 
family of Markov processes that is analogous to the Bessel-squared family. 
We show that such a family exists for each /3 and we establish some of the 
properties of these families. 

Recall that the transition probabilities of the unconditioned (1 + /3)- 
stable continuous state branching process, {Xt)t>Q, are characterized by the 
Laplace transforms, 

E^[exp(-0Xt)] =exp(-x0(t0'^ + l)-^/'^). 

Hence the transition probabilities of the conditioned process Y are charac- 
terized by the Laplace transforms, 

Ey[exp{-eYt)] = -Ey[eM-OXt)Xt] 

y 

y oO 

thus establishing (1.2). 

Recall also that if /3 = 1, then {Xt)t>o and (Yt)t>o are, up to a constant 
multiple, the Bessel-squared processes with dimensions and 4, respectively. 
The Bessel-squared process, {Zt)t>o, with dimension d (not necessarily inte- 
gral) is (up to constants) the Markov process characterized by the Laplace 
transforms, 

E^[exp(-0Zt)] = exp{-z9{t9 + ir^){t9 + ly^^^ 

This suggests that for < /3 < 1 and 5 > 0, there might be a semigroup, 
(Pt)t>o, such that 

(6.1) Pt exp{-9-){x) = Qx.^{-xe{te^ + i)-^'P){teP + ly^. 

We first verify that, for a fixed value of x, the right-hand side of (6.1) is 
the Laplace transform of a probability distribution (as a function of 9). We 
already know that 

eM-x0{t9^ + ly^^^) 
is the Laplace transform of a probability measure, so it suffices to show that 

{t9^ + 1)-"^ 

is also a Laplace transform of a probability distribution. Let {St)t>o be the 
/3-stable subordinator starting from 5o = normalized so that 

E[exp(-eSf)] =exp(-0'^i) 
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and let {Tt)t>o be the gamma subordinator starting from Tq = normalized 
so that for t>0, 

yt-l 

P{rtedy} = |^exp(-y)dy, 

and hence 

E[expi-eTt)] = {e + l)-'. 
Then if S and T are independent, 

E[eM-dStTs)] = neM-O^tTs)] = {tB" + 1)"^ 
We next verify that {Pt)t>o is a semigroup. Observe that 

= {te^ + l)-^Ps exp(-^(te^ + l)-^/^-){x) 

= exp(-xe((s + t)e'^ + + 1)^'' + 1)"^ 

= Ps+texp{-9-){x). 

It is clear that lim^o -ft exp(— 6'-)(x) = exp(— 0x) and so lim^o -fi/(a;) = 
f{x) for /gCo(M+). Also, lim^^a, Pj exp(-6'-)(y) = exp(-6'-)(x), and so 
limy_5.a; Ptf{y) = f{x) for / G Co(M+). The standard Feller construction gives 
that there is a strong Markov process, {Zt)t>o, with semigroup {Pt)t>o- 

This family of Markov processes shares many features of the Bessel- 
squared family. For example, it follows for a,b> 0, that 

Thus, the process {b~^^^ Zf,t)t>o is Markovian with the same transition prob- 
abilities as Z. Similarly, if Zq = 0, then the process (e~*/^Zgt)tgK is Marko- 
vian and stationary. 

Furthermore, if (Z^')f>o and {Z")t>o are two independent such processes 
with parameters 6' and 6" , then the process {Z^ + Z")t>o also belongs to the 
family and has parameter 5' + 6" . 

In a forthcoming paper, we will present a more thorough study of this 
family along the lines of [27, 28]. 
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